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We analyze the form of static charge susceptibility x(l) in underdoped cuprates near axial mo¬ 
menta (Q, 0) and (0, Q) at which short-range static charge order has been observed. We show that 
the momentum dependence of x(s) is anisotropic, and the correlation length in the longitudinal 
direction is larger than in the transverse direction. We show that correlation lengths in both di¬ 
rections decrease once the system evolves into a superconductor, as a result of the competition 
between superconductivity and charge order. These results are in agreement with resonant x-ray 
scattering data [R. Comin et al., Science 347, 1335 (2015)]. We also argue that density and current 
components of the charge order parameter are affected differently by superconductivity - the charge- 
density component is reduced less than the current component and hence extends deeper into the 
superconducting state. This gives rise to two distinct charge order transitions at zero temperature. 


Introduction.- Understanding charge order (CO) in 
high-T c cuprates and its interplay with superconductivity 
is essential for the understanding of the complex phase di¬ 
agram of these materials. An incommensurate charge or¬ 
der, accompanied by spin-order, was originally discovered 
in La-based cuprates^, but recently was also found to 
occur in Y- Bi-, and Hg- based materials^—, without an 
accompanying spin-order. A true long-range charge order 
has so far been observed only in a finite magnetic field^, 
but short-range static order (probably pinned by impu¬ 
rities) has been detected already in zero field. The CO 
has an axial momentum Q = Q y = (0, Q ) or Q x = (Q , 0) 
with Q ~ (0.2—0.3) x 2-7T, and can potentially be uni-axial 
(stripe), with only Q x or Q y within a given domain, or 
bi-axial (checkerboard) with Q x and Q y present in every 
domain. Recent STM and x-ray experiment s— 41 point 
towards the uni-axial order, at least at small values of the 
doping. The CO is often termed as charge-density-wave 
(CDW) to emphasize that it develops with incommensu¬ 
rate momenta, although its on-site component is sublead¬ 
ing to bond component because the measured form-factor 
for CO has predominantly a d-wave form£. 

The origin of the CO is still a subject of in¬ 
tense debate s 2 i 3 i 12 ~ 14 i 16 i 17 i 19 ~ 2i. Within one scenario, 
the charge order is induced by soft antiferromagnetic 
fluctuations^^ 14. This fits into the generic scenario that 
antiferromagnetism is the primary order parameter (i.e., 
the one whose fluctuations develop already at “high” en¬ 
ergies, comparable to the bandwidth), while CDW, its 
cousin pair-density-wave (PDW), and uniform d -wave su¬ 
perconductivity all develop as secondary orders induced 
by soft, low-energy magnetic fluctuations before the sys¬ 
tem becomes magnetically ordered. In another scenario, 
charge order is induced by lattice vibrations^; in this 
case lattice and electronic degrees of freedom should be 
taken into account on equal footing. And in yet an¬ 


other scenario, CO emerges in the process of the sys¬ 
tem transformation from a conventional metal to a Mott 
insulator—. If charge order reflects the crossover towards 
Mott physics, then the tendency towards localization of 
electronic states cannot be neglected even above optimal 
doping. 

A way to distinguish between these scenarios is to use 
the existing experimental data, particularly the ones for 
which the data analysis does not involve fitting parame¬ 
ters. Recent x-ray scattering data in underdoped YBCO, 
reported in Ref. eh can be used for this purpose. The 
data shows that the momentum structure of the charge 
susceptibility xifi) near Q x and Q y is anisotropic, and the 
longitudinal correlation length is larger than the trans¬ 
verse one, i.e. if \(g) at q = Q + q is approximated by a 
Lorentzian x -1 (9) ~ £ -2 + A 2 g| + A\q\, then Ay > A±. 

In this paper, we verify whether this condition is repro¬ 
duced within the two itinerant scenarios - the magnetic 
one and the phonon one. Exploring the Mott scenario is 
beyond the scope of this work. 

The charge susceptibility of itinerant electrons, x(q)> 
is generally related by a RPA-type formula to the static 
particle-hole polarization bubble between low-energy 
fermions separated by q (“hot” fermions). The differ¬ 
ence between the two scenarios is that in the magnetic 
one the interaction is peaked at momentum K = ( 7 r, 7 r) 
and connects fermions from two different hot regions. 
Specifically, magnetic interaction moves the center of 
mass momentum of a pair of hot fermions from ko to 
= ko + K. Then, one needs to apply spin-fluctuation 
mediated scattering twice to move fermions back to the 
same hot region (see Fig. [TJ, and, as a consequence, 
x -1 (?) oc l-U 2 n fco (q)n fc)r (q)> where II fco (q) and n fcir (q) 
are polarization operators made out of hot fermions with 
relative momentum q and center of mass momentum near 
ko and k n , respectively (see Fig. Or). In the phonon see- 
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nario, the interaction acts independently within each hot 
region, and x _1 (g) oc 1 - £7[II feo (q) + II fcjr (q)]. In both 
cases, the momentum dependence of x(q) is solely deter¬ 
mined by the polarization bubbles and does not depend 
on the strength of the interaction U or U. 

The polarization bubbles IIfc 0 (q) and IIfc 7r (q) depend 
on the Fermi surface geometry in the vicinity of the hot 
spots and also on the choice of the upper cutoff in the 
momentum deviations from the hot spots. The presence 
of the cutoff, A, reflects the fact that the interaction in 
the charge channel can be approximated by a constant 
(U or U) only in a finite range around a hot spot, out¬ 
side of which it drops rapidly. In particular, within spin- 
fluctuation scenario, A depends on the distance to the 
magnetic QCP it tends to a constant at the QCP and 
scales as inverse magnetic correlation length when 
£ s drops below a certain valu o 24 i 27 . 

We report the results of analytic computations of 
IIfc 0 ( q ) and 11*^ (q) using a hard cutoff (to be defined be¬ 
low). In the magnetic scenario, the longitudinal charge 
correlation length turns out to be larger than the trans¬ 
verse one, in agreement with the data^i. For the phonon 
scenario, the result is the opposite - the transverse cor¬ 
relation length is larger. Taken at face value, this obser¬ 
vation selects magnetic mechanism of CO formation over 
the phonon one, although it indeed does not preclude 
phonon scattering as a subleading mechanism of CO. 

We also consider how charge susceptibility x{ ( i) gets 
modified once the system becomes a d- wave supercon¬ 
ductor. We find that the key effect of superconductivity 
is the reduction of IIfc 0 (Q) and IIfc 7r (Q) due to competi¬ 
tion between CO and superconducting order parameters. 
As a result, both longitudinal and transverse CO corre¬ 
lation lengths get smaller in the superconducting state. 
A more subtle result is that the reduction of the bub¬ 
ble is different for density and current components of 
CO. The two are symmetric and antisymmetric combi¬ 
nations of the incommensurate charge order parameters 
Afc = c k+Q/ 2 , a S <*P c k-Q/ 2 ,P witl1 k = ±k o- The current 
component changes sign under time-reversal and once it 
develops along with the density component, the CO spon¬ 
taneously breaks time-reversal symmetry^^. In the nor¬ 
mal state, density and current susceptibilities are equal 
[both are x(<?)], as long as k 0 and — k 0 are well separated 
such that one can neglect bilinear coupling between A^ o 
and A5 fco . In the superconducting state, X -1 (Q) for the 
current component of CO gets shifted by A 2 sc /{TKvf) 
due to negative feedback from long-range superconduct¬ 
ing order, A sc . For the density component of CO, such 
term cancels out and the shift is much smaller, of order 
A 2 C /(A vf ) 2 x log(A/T). As a result, the density com¬ 
ponent of CO is much less affected by superconductiv¬ 
ity and should persist deeper into the superconducting 
state. This is in agreement with the x-ray data, which 
found that the measured charge density fluctuations per¬ 
sist down to lowest temperatures^, deep into the su¬ 
perconducting state, where they also likely get pinned 



FIG. 1. The Fermi surface (black) with the occupied states 
shown in green. The location of the hot spots for CO (the 
points on the Fermi surface separated by Q y = (0, Q) or Q x = 
(Q,0)) are marked as 1, 2, —1, —2,... with the corresponding 
directions of Fermi velocities, vf■ The notations ±&o and 
±fc,r are for the center of mass momentum of charge order 
parameter A® = c\ +Q/2 J a pc k . Q/2 ^. 



FIG. 2. The Feynman diagrams for (a) IIfc 0 (O) and Ylk„(Q) 
in the normal state, and, (b) normal and anomalous contri¬ 
butions to II feo (Q) in the superconducting state, (c) Four- 
point diagrams contributing to /3ii and /J 12 in the free-energy 
(Eqn[2| in the normal state. 


by quenched disordei—. Likewise, at zero temperature, 
the CO emerging from pre-existing superconducting state 
should initially have no current component. 

Polarization operators in the normal state.- For the 
computation of polarization operators we used the Fermi 
surface shown in Fig. [T] For definiteness we set Q = 
Qy = (0,Q). The results for Q = Q x are identical by a 
7t/2 rotation. There are four “hot” regions in the Bril- 
louin zone. For the two regions with ko = (n — <3/2,0) 
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and — ko, the Fermi velocity of the two fermions sepa¬ 
rated by Q are almost anti-parallel, while for the other 
two, with k^- = {—Q/ 2,7t) and — k^, the velocities are 
nearly parallel. In the notations in Fig. [lj this implies 
v y v x ( v v /v x « 13.6 in BSCCO, see RefTl30t). 

We impose hard cutoff by requiring that momentum 
of each fermion in the particle-hole bubble IIfc 0 (g) and 
II k„{q) differs from the corresponding hot spot by no 
more than A by amplitude. In the hot region with 
center of mass momentum of a pair near we define 
ki ,2 = k 0 ± (Q y + q)/2 + k, and the condition reads 
|k ± q/2| < A. The analogous condition holds in the 
hot region with center of mass momentum of a pair near 
k n . We assume that A is small compared with inverse 
lattice spacing, in which case we can expand the disper¬ 
sion of a hot fermion to linear order in deviation from 
a hot spot. Under these condition, we obtained analyt¬ 
ical expressions for the polarization bubbles rU 0 (g) and 
II k„{q) to leading order in q = q — Q y . The calculation 
is lengthy but straightforward^, and the result is 

iWq) = Av x (2 A log _ \q y \ sin" 1 ^ - \q x \log , 

Ufc„ (q) = Av y (2Alog ^ - \q y \ sin" 1 ^ - \q x \ log 

(1) 

where v = ^Jv 2 + v 2 and A = l/(n 2 v x v y ). It is interest¬ 
ing to note that the singular terms proportional to \q x , y \ 
are independent of the cutoff, A. 

The case of hard cutoff is somewhat special because the 
expansion of polarization bubbles in q is non-analytic and 
holds in powers of \q x \ and \q y \ rather than in q 2 and q 2 . 
We verified that the quadratic dependence emerges imme¬ 
diately once we soften the cutoff. Still, even for strictly 
hard cutoff one can analyze the anisotropy of the inverse 
charge susceptibility X” 1 ( ( ?) by comparing the prefactors 
for longitudinal (|g y |) and transverse (1^1) momentum 
dependencies. By continuity, the anisotropy should sur¬ 
vive upon softening of the cutoff. 

To analyze the anisotropy, we use v y v x and 
expand IIfc 0 and 11^ in small v x /v y limits. In the 
magnetic scenario we find X -1 (9) Co + C y \q y \ + 
C x \q x \ + ..., where C 0 = 1/U 2 - 4A 2 A 2 i> 2 log 2v y /v x , 
Cy = 2A 2 Av 2 (log 2v y / v x +7r/2),C x = 2A 2 Av 2 logv y /v x 
and the ellipses denote higher order terms in q x , y - We 
remind that q = Q y + q, hence q y is longitudinal com¬ 
ponent. Taking the ratio C y /C x we immediately see 
that C y /C x = 1 + [(7r/2 + log2)/log (v y /v x )\ > 1, i.e., 
the effective correlation length £|| = C y /C 0 is larger 
than = C x /Cq. This is consistent with the data^. 
For v y /v x « 13.6, we obtained, without expanding, 
£||/£_l = 1-87, which is reasonably close to the experi¬ 
mental ratio of around 1.5. 


In the phonon scenario we obtain \ 1 (o) °c Co + 
C y \q y \ + C x \q x \ + ..., where now, to logarithmic accu¬ 
racy, C v /C x = (ir/2 + l)/log (v y /v x ). Then C y /C x is 
smaller than one, at least when v y /v x is large enough. 
For v y /v x « 13.6, we obtained, without expanding, 
?||/a = 0.98. 

We also computed nfc 0 (q) and 11/^ (q) numerically for 
a specific Lorentzian cutoff, which we imposed by insert¬ 
ing into the integrands for the bubbles an additional fac¬ 
tor, A 2 /(A 2 + (k + q/2) 2 ) x A 2 /(A 2 + (k - q/2) 2 ), but 
not restricting integration over momentum. One can im¬ 
mediately make sure that in this case the expansion in q 
holds in powers of g 2 . We again find that in a magnetic 
scenario £||/£_l > 1- However this ratio is much larger. 
From this perspective, the hard cutoff gives better agree¬ 
ment with the data. 

Superconducting state.- The polarization operators 
and n^g) in the superconducting state are ob¬ 
tained in a conventional way, by combining bubbles made 
out of normal and anomalous fermionic Green’s functions 
(Fig. [2jo) . The full expressions are more involved and we 
didn’t obtain analytical formulas even for hard cutoff. 
In general, both n/ Co (Q y )/Hfc 7r (Q y ) and the momentum- 
dependent terms in the polarization operators evolve 
with the superconducting gap A sc . The effect, however, 
is stronger for y(Q y ) rather than for g-dependent terms 
because for x(Q y ) superconductivity-induced shift has to 
be compared with the initially small value of the mass of 
the charge susceptibility at Q y . We therefore focus on 
the renormalization of the polarization operators right 
at q = Q y . The calculations, which we describe in more 
detail below, show expected trends - superconductivity 
competes with CO, and once long-range superconduct¬ 
ing order develops, it tends to delay the appearance of 
CO. This effect is very typical for competing orders and 
has been recently discussed in detail for Fe-pnictides^. 
Because X -'(Qy) increases, both longitudinal and trans¬ 
verse charge correlation lengths go down. The data*ii 
show the same trend. 

A more subtle issue is the magnitude of 
superconductivity-induced shift. Near T c (i.e., 
for relatively small A sc ), the shift originates from 
f3ij\A sc \ 2 A®. (Aj?)* terms in the Free energy, where, we 

remind, A® is fluctuating CO field (not the condensate), 
and hi is either ±fco or ±fc T . Since the superconducting 
pair has zero total momentum, it couples to fermions 
located only within one corner of the Brillouin zone; 
hence coupling terms for ±fco and ±fc 7r can be considered 
separately. At the same time, superconductivity pairs 
fermions with opposite momenta, hence both |A®| 2 and 
A^(A^ fe )* couple to |A SC | 2 . 

For definiteness, let’s set, as before, Q = Qy and focus 
on the region where center of mass momentum of charge 
order parameter is ±fcg. The Free energy is 


F = X - 1 (Q y ){\A^\ 2 + \A Q _l o 


+ |A SC | 2 [/3 n (|Agf+ |A?^ 


I 2 ) +P 12 (a£*(A %)* + {A%)*A% 


( 2 ) 
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(We neglected spatial fluctuations of CO and terms un¬ 
related to our purposes.) This Free energy is easily di¬ 
agonalized by introducing A d = (A^/ + A®£ o )/\/2 and 
A c = (A^ — A® v ko )/y/2. In terms of these variables 

F= |A d | 2 (x-^Qy) + 2|A sc | 2 (/3n + ft 12 )) 

+ |A C | 2 (x-^Qy) + 2|A sc | 2 (/3 n - ftu)) . (3) 

We see that the shift of X~ 1 (Qy) due to superconduc¬ 
tivity is generally different for symmetric density and 
antisymmetric current components of CO (Ad and A c , 
respectively). 

The couplings fa and fa can be evaluated either by 
using Hubbard-Stratonovich (HS) formalism, or by ex¬ 
panding particle-hole bubbles to order A 2 C . In the HS 
formalism, these two terms are given by square diagrams 
made out of four fermionic Green’s functions in the nor¬ 
mal state (shown in Fig. [2t). 

We computed fa and ft 12 f° r a model with hard cutoff 
and found that the dominant piece in each is a cutoff- 
independent term, which scales as 1/T. At the lowest 
T, the expansion in powers of A sc does not hold, and 
1/T divergence is cut by 1/|A SC |. Upon a more careful 
look, we found that universal (i.e., cutoff-independent) 
1/T terms in ftu and ftu come with exactly opposite co¬ 
efficients, i.e., 1/T terms cancel out in ftu + ft 22 . This 
cancellation has not been noticed before. From math¬ 
ematical perspective, this cancellation is similar to the 
cancellation between vertex and self-energy corrections 
to uniform density-density correlator. In the latter case, 
however, the cancellation is exact and it enforces Ward 
identity associated with particle number conservation. In 
our case, the cancellation is not required by symmetry 
and is not exact - only the leading 1/T terms cancel out 
in ftu + fti 2 . The subleading terms do not cancel. These 
subleading terms are, however, much smaller, and the 
correction to X^iQy) from the superconducting order is 
given by a|A sc | 2 /A x log(A/T) (a > 0). The outcome 
is that the mass of the density component of charge or¬ 
der parameter (the one measured by x-ray) goes up in 
the presence of superconductivity (and the correlation 
length, which scales as inverse mass, goes down); however 
this effect is small. A small coupling between competing 
charge density and superconducting orders also implies 
that the order which appears first does not rapidly de¬ 
stroy the other one, hence CDW and superconductivity 
co-exist over a sizable range of dopings. 

At the same time, for the current CO component A|?, 
A 2 C /T contributions in ftu and ftu add up, i.e., for this 
component superconductivity has stronger negative im¬ 
pact. In particular, when CO emerges inside the super¬ 
conducting dome, only its density component becomes 
non-zero when -x -1 (Q) = | A sc \ 2 {ftu +^ 22 ) ~ |A SC | 2 /A. 
Current component emerges at smaller dopings, when 
(and if) the condition -~X~ 1 (Q) = |A sc | 2 (/3n - ft 22 ) ~ 
|A sc | 2 /(max(T, |A SC |)) is satisfied, and co-exists with su¬ 
perconductivity in a narrower range. We show this in 
Fig. [3| Since the current component of CO is respon- 



FIG. 3. The phase diagram in variables T and x (adapted 
from Ref. H~5h . with more details near T = 0. N and TRSB 
stand for nematic and time-reversal symmetry breaking, re¬ 
spectively, and CDW is the phase with broken translational 
symmetry (in a clean system). 


sible for the breaking of time-reversal symmetry (TRS), 
its absence over some doping range where the density 
component of CO is present implies that in this range 
CO cannot be the source of TRS breaking. Alternatively 
speaking, if TRS breaking is caused by charge order, the 
end point of TRS breaking transition should end up at 
T = 0 at a smaller doping than the onset of CO. It would 
be interesting to test this in Kerr and elastic neutron scat¬ 
tering measurements in the superconducting stated—. 

Summary.- In this paper we considered three aspects 
associated with uni-axial charge order in the underdoped 
cuprates. First, we analyzed the anisotropy of the charge 
order correlation length in the normal state, detected 
in recent x-ray measurements. Our goal was to investi¬ 
gate whether these data allow one to distinguish between 
magnetic and phonon-based mechanisms of CO forma¬ 
tion. We argued that the magnetic scenario yields results 
consistent with the data in Ref. E0. Second, we argued 
that both longitudinal and transverse charge correlation 
lengths decrease in the presence of a true superconduct¬ 
ing order, primarily because this order increases the mass 
of the charge order propagator. Finally, we found that 
the mass increase is different for the density and the cur¬ 
rent components of CO (symmetric and antisymmetric 
components with respect to the flip of a center of mass 
momenta of fermions which form CO). The mass increase 
is strong for the current component and is parametrically 
weaker for the density component. As a result, under the 
umbrella of superconductivity, the density component of 
CO exists in a wider range of doping compared to the 
current component. Since the current component of CO 
is responsible for TRS breaking, we propose to use this 
fact to test whether TRS breaking is associated with in¬ 
commensurate charge order. 
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Supplementary Material 


I. EVALUATION OF THE PARTICLE-HOLE BUBBLES IN THE NORMAL STATE 

In this section we evaluate the particle-hole bubble IIfc 0 and associated with the incommensurate charge order 
with momentum around, say, Q = Q y . The results obtained below can be trivially extended to the case for Q = Q x 
by rotational symmetry. 

By simple dimensional analysis the particle-hole bubble ~ f diocPk GG , where G is the Fermionic Green’s function 
in the normal state, scales as the upper momentum cutoff A. For the spin fluctuation scenario, the cutoff scheme 
in momentum is natural since the structure of the interaction function dictates that only fermions in the hot region 
are relevant. For a small antiferromagnetic correlation length £ s , this cutoff around corresponding hot spots is given 
by roughly A ~ £“ * 1 . On the other hand, for the phonon scenario, it is not clear how to impose the momentum 
cutoff through microscopic mechanism. However, such a cutoff near hot regions is still phenomenologically required - 
otherwise the local minima of n(Q) will be located at Q’s connecting the nested portions of the Fermi surface near 
the antinodes, rather than the experimentally-detected^ Q = Q y . To ensure an unbiased comparison between the 
two scenarios, we confine the typical momentum deviation from hot spots to be of order A, without specifying the 
microscopic origin of A. We adopted both “hard” and “soft” cutoff schemes, and found that a hard cutoff yields a 
better match with the experimental data. 


A. Particle-hole bubbles with a hard cutoff 


We now evaluate the bubble near hot spots 1 and 2 with a hard momentum cutoff. 

d 2 k 1 1 

'|k±q/2|<A 


Iffeo (Q + q) — —T ^2 


(2 tt) 2 


- ei(k + q/2) iio m - e 2 (k - q/2)' 


( 4 ) 


We first evaluate this integral for q = q x . Performing the summation in frequency we have, 


nfe 0 (Qx) 


' |k±q/2|<A 


d 2 k /(ei(fc + q x / 2)) - /(e 2 (fc - q x / 2)) 
4-7T 2 e 1 (k + q x /2)-e 2 (k-q x /2) 


J_ f d 2 k 

47r 2 J c \2vyky + v x q x \' 


( 5 ) 


where /(e) is the Fermi function. In the last step we have approximated /(e) by the step function, and we have 
used the linear dispersion ei j2 (fc) = v x k x ± v y k y . The integration region C is given by {ei (k + q x /2)e 2 (k — q x / 2) < 
0}n{|k±q/2| < A}, which we depict in Fig.[U)b) for the case v x \q x \ <C A. For comparison we also show the integration 
region for q = 0 case in Fig. 0)(a). 

We can then divide the integration region into more regular shapes and evaluate them separately. Before we do so, 
it is convenient to introduce rescaled variables y = 2v y k y /(v x q x ), x = 2k x /q x , and Q = 2v y K/(v x \q x \): 


n fco ( q x ) 


\q x \ f dx dy 

16n 2 v y J c , \y+l\' 


( 6 ) 


where C is now given by {|x| < |y + l|}n{|x| < ^yjQ 2 — y 2 — 1}- We take the limit Iq^l <C A, and hence Q> 1. By 
dimensional analysis, the integral from Eq. JG]) in leading order scales as Q, which, when combined with the prefactor 
|g x |/(167r 2 u y ) gives a contribution ~ A. This is nothing but n/ Co (Q). To obtain the leading order correction for 
q = q x , one needs to keep contributions up to next leading order in Q. 
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FIG. 4. The integration regions for particle-hole bubble near hot spots 1 and 2 with (a) total momentum right at Q, (b) 
momentum deviation q = q x , and (c) momentum deviation q = q y . The circles centered around ±q/2 with radius A are 
imposed by the hard upper momentum cutoff |k ± q/2| < A, and the straight lines with slope izv x /v y are FS’s in the vicinity 
of hot spots 1 and 2 . 


Using Fig. QJb), we can divide the integration region C' as IIfc 0 (g x ) = \q x \/(\d>TT 2 v y ) (^f c , + f c , + f c , 
\q x \/(l&Tt 2 v v ){h +I 2 + h), where 


h 

h 

h 


dx dy 

\y + 1 | 

dx dy 

I y + M 

dx dy 

\y +1| 


r Qv x /v-2v 2 /v 2 
' -Qv x /v 

r\JQ 2 - v ll v l 


dy 


dy 


r\y+l\ 1 

/ dx- -— 

-|y+l| |Z/ + 1| 

•v x ^Q 2 -y 2 /Vy~l 1 

dx- 


'Qv x /v—2v^/v 2 J-V x y/Q 2 -y 2 /vy + 1 12/ ~F 11 

r—QVx/v 


dy 


r:\jQ 2 -y 2 hy-l 

dx 7 


-V Q 2 -»ihi J-vx^/o^/vy+i I y + M 


(7) 


The integration limits for I\ comes from the condition \y + 1| < ( Vx/v v )\ / Q 2 — y 2 — 1 and has been expanded up to 
next leading order in Q. 

The evaluation of I\ is trivial, which, to next leading order in large Q , gives 

h = 4 Qv x /v - 4 vl/v 2 . ( 8 ) 

For I 2 + I 3 we have 



In the second line above we have kept terms only up to next leading order in Q. Note, however, that the above result 
is not an approximation in the small v x /v y limit but rather for an arbitrary v x /v y . 

Combining the results for Ii. 2 , 3 -, and using II k 0 (qx) = |< 7 a;|/(167T 2 (Ti + J 2 + I 3 ), we find in original variables 


n fco (qx) 


A v + Vy \q x \ v 

&s log — 


( 10 ) 


The first term of Eq. COD corresponds to IIfc 0 (q x = 0) = IIfe 0 (Q). This term has been evaluated^ before, with a slightly 
different cutoff scheme. In the limit v y v xi we have IU o (0) = A/(2tt 2 v v ) log (vy/vx), and this is fully consistent with 
previous result. 












































Next we evaluate the same bubble with momentum deviation q = q y from the CDW momentum Q. Similar to Eq. 
© we start from 


Hfco(9y) — / 

J li 


d 2 k f(e\(k + q y j 2 )) - f(e 2 (k - q y / 2 )) 


|k± q /2|<A 4tt 2 ei (k + q y /2) - e 2 (fc - q y /2) 
d 2 k 


( 11 ) 


47T 2 J d |2 v y ky\' 

The integration region D is given by {ei(/c + q y /2)e2{k — q y / 2) < 0} fl {|k ± q/2| < A}, which we depict in Fig. dfc) 
for the case v y \q y \ <C A. Rescaling y = 2 k y /(q y ), x = 2v x k x /(v y q y ), and Q = 2v x A/{v y \q y \) 1, we rewrite Eq. (fill) 
as 


Hfco(9y) — 


\%\ 


167r 2 z;. 


x J D' 


dx dy 

\y\ 


( 12 ) 


For simplicity of presentation, we drop the bars in the rescaled variables below. From Fig. [3Jc) we see that the 
integration region D' can be simply expressed as \x + 1| < |y| < ( v y /v x )y/Q 2 — x 2 — 1, where Q > 1. We then 
integrate over y first and obtain 


Hfco (%) — 


Wvl [ {VjV)Q ~ 2VI/V2 dxlQg ( Vy_/V X )VQ 1 


8tT 2 V- 


* •> ~(Vy/v)Q 


X + l 


(13) 


where the integration limits on x is obtained by expanding the solution of \x + 1 | = ( v y /v x )\/Q 2 — x 2 — 1 up to next 
leading order in Q. Evaluating Eq. (1131) in the large Q limit, we obtain 


Hfco (%) — 


m\ 


8tt 2 v x 


[ (Vv/V)Q Jxlo E (^)V^ [ iVy/V)Q 

'-(Vy/v)Q l^l 


r(v y /v)Q dx /-0 


dx 


{v y /v)Q-2vllv 2 
( v y/ y )Q dx 


dx log 


{v y /v x )yjQ 2 - x 2 


'—( V y/V)Q \ X \ V U J-{Vy/v)Q Q 2 — X 2 


(14) 


where the first term in the bracket is of order Q , which corresponds to IIfc 0 (Q), and the last four terms are of order 
one. The second Evaluating the integrals, we obtain 


TT / \ A V + Vy 

I W*) = SS - 1 ° g — 


kv I . 


-i u y 

— 5 —sin —. 

47 ^ 1 !,, v 


Combining Eqs. (111)1151) . we get the full expression 

_ , _ , A v + v. 

n fco (Q + q) = 7 ^— log 


2tt 2 i 


y i 9 x| , v \q y \ . _ 1 v y 
log- , " sin 


4tT 2 V v 


47r 2 i;„ 


(15) 


(16) 


Finally, we move to the evaluation of 11^, defined as the particle-hole bubble for fermions near hot spots 3 and 4. 
Similar to Eq. ®, we have 


(Q + q) - —T X] 


d 2 k 


1 


1 


|k±q/2|<A ( 27r ) 2 iUm - e 3 (k + q/2) iu m - e 4 (k - q/2)' 


(17) 


We remind that e 3i 4 = —v y k x =F v x k y , while ei j2 = v x k x ± v y k y . Using the fact that the integrations over k x and k y 
are symmetric, we then see that IIfc 7r (Q + q) and IIfc 0 (Q + q) are simply related by v x -O- v y , therefore 


(Q + q) 


A v + v x \q x \ v \q y \ . _ x v x 

-is;"" - 


(18) 


As a quick check, in the limit v x <C v y , the q independent term becomes 11^ (Q) = A/(2n 2 v y ). This is consistent 
with results previously obtained^. Eqs. (I16I18D are the two main results of this Subsection. 

When v y v x we can expand IIfc 0 and II^ in small v x /v y , 


n fco (q) « A 
(q) « A 


7T V 

2Av x log —- - \q y \-^v x - \q x \v x log -A 
Vx ^ V x 

Ip 

2A.V X lO'yl'^a; \Qx | T 

ZVy 


( 19 ) 
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B. Particle-hole bubbles with a soft cutoff 


Using a hard cutoff, we found that the leading order q dependence of the polarization bubbles has a non-analytical 
form ~ |g|. This |g| dependence gets “softened” to a q 2 one, if we use a soft upper cutoff instead of a hard one. For 
example, one can replace the integration limit (k±q/2 ) 2 < A 2 with a multiplier A 2Ar /(k + q/2) w +A Ar )/(k — q/2) w + 
A N ) and integrate over infinite range. This multiplier is fully analytical, hence the leading order q dependence should 
be ~ q 2 . For large N, this multiplier almost looks like a step function, which restores the hard cutoff. As a result, 
the coefficient of q 2 term increases, which diverges at N —> oo and gives a ~ |g| behavior in the hard cutoff case. We 
have verified that the ratio between coefficients for q 2 and q 2 for a large but finite N connects smoothly to the ratio 
between \q x \ and \q y \ coefficients obtained previously using a hard cutoff. 

For completeness, in this subsection we perform a calculation of the bubbles for N = 2, which is a Lorentzian cutoff. 
Namely, we define 


Hfco (Q 

iMQ- 


<0 = -r£ 

m ' 

q) = - r E 


cPk 


1 


1 


A 4 


(27r ) 2 iu m - ei(k + q/2) iu m - e 2 (k - q/2) [(k + q/2) 2 + A 2 ][(k - q/2) 2 + A 2 ] 

d 2 k _1_1_A 4 _ 

(27r ) 2 iu m - e 3 (k + q/2) iuj m - e 4 (k - q/2) [(k + q/2) 2 + A 2 ] [(k - q/2) 2 + A 2 ]' 


( 20 ) 

( 21 ) 


Using £ 1 , 2 (k) = v x k x ± v y k y , we can redefine k = (k x , k v ) = ( k x + v y /(2v x q y ), k y + v x /(2v y q x )), such that in Eq. 
(tZHD we have ei (k + q/2 ) = ei(fc) and e 2 (fc — q/2) = e 2 (fc). This way all the q dependences are absorbed into the 
Lorentzian multipliers: 


IIfc 0 (Q + q) — — 


d 2 k /(ei(fc)) - f(e 2 {k)) 


A 2 


47 r 2 ei(k)-e 2 {k) 


[kx Vydy/{2v x ) q x /2] 2 -f- [ k y v x q x /i/^Vy) T <?y/2 ] 2 U A 2 

A 2 


[k x v y q y /(2v x ) qx/2)/ 2 T [ k y v x q x /(2v y ) q y / 2 ] 2 T A 2 


( 22 ) 


where, as before, we have summed over frequency. To obtain the quadratic coefficients in q' s, we can expand the 
Lorentian multipliers in q x and q y s. The result will then have the form II feo (Q + q) = II fco (Q, r) +D x (r)q 2 +D y (r)q 2 , 
where we keep the r = v y /v x dependence explicit, and 


D x 


d 2 k /(ei(fc)) — f(e 2 (k)) d 2 


8tt 2 e L (k) — e 2 (k) 


A 2 


dQ x ( [k x v y q y /(2v x ) T q x / 2 ] 2 + [ k y v x q x /(2v y ) -t- q y /2[ 2 

A 2 


A 2 


[k x v y q y /{2v x ) q x /2[ 2 T [ k y v x q x /(2v y ) q y /2\ 2 A 2 


9=0 


D v 


d 2 k /(ei(fc)) — f(e 2 (k)) d 2 
8 tt 2 ei(fc) — e 2 (k) dq 2 


A 2 


[k x v y q y /(2v x ) T q x / 2 ] 2 T [k y v x q x /(2v y ) T q y / 2]“ T A~ 


A 2 


[k x ^yQy/(/^Vx) q x /2i\ 2 -f* [k y v x q x j(2v y ) q y /2[ 2 T A 2 


(23) 


9=0 


Like we did before, 11*^ and IIfc 0 are related by k x O k y , hence h^ (Q + q) = IIfc 0 (Q, 1/r) + D x (l/r)q 2 + D y (l/r)q/. 
In the magnetic scenario, the spatial anisotropy of CDW correlation length £||/£_l is given by the ratio of q 2 coefficients 
in IIfc 0 (g)!!^ (q), which means 


C||_ = n 0 (Q, r)D y (l/r) + II 0 (Q, l/r)D y (r) 
n 0 (Q, r)D x (l/r) + n 0 (Q, l/r)D x {r )' 


(24) 


We used the experimental ratio r = 13.6 and evaluated all the coefficients numerically, and we obtained £||/£_l = 27.3. 
This ratio, however, is much larger than the experimental value ~ 1.5. 


II. COUPLING BETWEEN SUPERCONDUCTIVITY AND CHARGE ORDER 


We next compute the effect of feedback from the superconducting (SC) order parameter A sc to the density and 
current component of the charge order (CO) parameter. We use two approaches. First, we compute the coupling 
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between CO and SC perturbatively by evaluating the square diagrams (see Fig.2(c) in the main text). Second, we 
compute the particle-hole bubble in Narnbu space, which naturally encapsulates the contribution from the SC order 
parameter A sc into the normal and anomalous Green’s functions, and at the end expand the result in terms of the 
SC order parameter. We show in Section Hi A1 using both methods, which are equivalent by nature, that the leading 
contributions cancel for density component, while they add up for current component. This leads to a much weaker 
“residual” feedback effect from the SC order parameter on the density component compared to the current component, 
which we compute in Section III Bl 


A. Cancellation of the leading order coupling between density component of charge order and 

superconductivity 


We start by doing perturbation theory. At lowest order the coupling between A sc and A’s are given by square 
diagrams shown in Fig. 2(c) in the main text. Note that there are two types of such diagrams. In the first type, only 
the modulus of the CO parameter with a given center of mass momentum, |Aj?| 2 , is involved. In the second type, 
however, the CO parameters involved are of the form A^(A® fc )*. Therefore they correspond to the coefficients /3u 
and /3i2 in Eq. (3) of the main text respectively. 

Like the bubbles, the square diagrams should also be evaluated for A’s with center of mass momenta at both 
±fco and +k n . However, the ratio between Aj? and A® is fixed since they are coupled by the antiferromagnetic 
interaction^. Moreover, by a simple transformation, one can show that the integrals for the square diagrams of the 
same type involving fermions near k 0 and k n are identical. Therefore, the contribution to the Free energy from these 
square diagrams involving Aj.^ just differ from those involving A±k 0 by a constant. For simplicity we will only 
consider the diagrams with A±k 01 and the extra constant can be easily absorbed into a rescaling of A sc . Therefore, 
we have 


P ii 
Pl2 


- 2T Y, I ^G 1 ( Wm ,A : )G 2 ( W7n ,fc)G 1 ( Wm ,fc)G_ 1 (-a; m ,-fc) 

m ^ 71 

v f rpk 

~ t ^2j ^2Gi(u m ,fc)G_i(-w m ,-fc)G 2 (u„„fc)G_ 2 (-w m ,-fc), 


(25) 


where Gi(w m , k) = 1 /[iu> m ~ £i(fc)]- Both integrals have been evaluated in Refs. [|] and@ with linear dispersion ei i2 
and infinite limits on momentum integration. However, an overall sign of /?i 2 was overlooked. The revised results are 
of the form, 


Pu 


P 12 


1 

32v x v y T 


1 + 0 



1 

32 v x v y T 


1 + 0 



(26) 


Plugging these into Eq. (4) of the main text, we immediately see that there is a cancellation of leading order (~ 1/T) 
contributions in the bi-quadratic coupling between the density component A^ and H>, while the bi-quadratic coupling 
between the current component A c remains ~ 1/T. 

The same cancellation between leading order contributions can also be derived within the Nambu formalism in the 
SC state. For simplicity we focus on the region of hot spots 1,2,-1,-2 only. We introduce the Nambu spinors as, e.g., 

+\(k) = (cj f (fc) c_ u (-fc)) , +\{k) = (c\ t {k) c_ 24 .(-fc)) . (27) 


The SC gap A sc can then be naturally included into the fermionic dispersion as 


'Ho + 'Hsc = '&\{k)[ei{k)T z + A sc t x ] +i(k) + +\{k)[e 2 {k)T z + A sc t x ]+ 2 (k), 


(28) 


where t x and t z are Pauli matrices in Nambu space, and for simplicity we have taken A sc to be real. From the above, 
we can write down the Green’s functions as, 


k) 

G-(w m , k) 
E?(k) 


_( G?(w m ,fc) 

y-F^^k) 

iu) m + ei(k) 

< + E?(ky 

=€ i(k) + |A SC | 2 . 


-F9(w m ,k) 

-G?(-w ro ,fc) 

F°(u m ,k) = 



‘—^SC 

< + Ef(k) 


(29) 


, where 


(30) 

(31) 
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We now couple the fermions T’s to charge order parameters A^ and A c . Using the fact that Ad( c ) is the even (odd) 
combination of the CO parameters a2^ 0 , we have 

Ha = A d (ty{(k)T z 'f) 2 (k)) + A c (^'{(fc)To^ , 2(fc)), (32) 

where to is an identity matrix. We evaluate the bubbles for A<j iC , and the coupling between A<j, c and A sc can be 
obtained by expanding these bubbles in powers of A sc . 

For Ad, we find 


IT 


fco 


(q.)=- t E 

m ' 

= 2 r E/ 


Tr[C/i (a; m i k)-T z ■ £ 2 (w m ,fc) ■ T z ) 


f drk 

(27r) 2 

d 2 k — ei(k)e 2 (k) 


(2 4 2 


■Ef(k)][^+Ei(k)Y 


(33) 


At this step if we expand in the integrand with respect to A sc to quadratic order, the coefficients we get are precisely 
f3n + Pi 2 - As we have seen, both /3’s scale as ~ 1/T and cancel out. Here we do something slightly different to show 
the same cancellation. Namely, we take the T —> 0 limit and replace the Matsubara sum with an integral, then 


nt(Qy) = 


du n 


d 2 k u 2 m - ei(k)e 2 {k) + |A S 


(27r) 2 [o& 


El{k)}[^ra 
d 2 k S 1 (fc)U 2 (fc)-ei(fc)e 2 (fc) 


-Eim 

|A SC | 2 


J (2 t)2 E 1 (k)E 2 (k)[Ei(k) + E 2 (k)] ' 

To compute the effect of A sc , we subtract from the above integral its A sc -independent part, namely 2Hfc 0 (Q). 

<sn&(Q„) = 


f d 2 k 

1" Ei{k)E 2 {k) — ei(k)e 2 (k) + |A SC | 2 

\ei(k)e 2 (k)\ - ei (k)e 2 (k) 

j (27r) 2 

E4k)E 2 {k)[E4k) + E 2 (k)} 

| ei (fc)e 2 (/c) | (| ei (/c) | +|e 2 (fc)|)_ 


(34) 


(35) 


To obtain leading order contribution, which we know comes from the infrared, we can evaluate Eq. (1351) with linear 
dispersion ei,e 2 and infinite integration limits. A little trick here is to transform the integration variables from 
k = (k x , k y ) to (ei, e 2 ). The Jacobian of this transformation is l/(2v x v y ). It is then clear that the eie 2 terms in both 
numerators are odd and hence can be eliminated. We regroup into two terms, 


sni(Q y ) = 


deide 2 


8tt 2 


V X V y 


\A S 


|A S 


8 V X Vy 


E\E 2 (Ei + E 2 ) 
00 d£ide 2 


1 


E\ + E 2 

1 


M + M 


71-2 \Jl + 4 + e 2 (\/l + 4 + \/l + e 2 ) 

d£\d£ 2 I 1 1 


7 T 2 l 1 -)- £ 2 + ^1 + 4 M + \£2 


(36) 


where we have defined £-\, 2 = ei >2 /| A sc |. Naively, this means the leading order renormalization to the bubble for charge 
order is of order ~ |A SC |. This |A SC | scaling is nothing but the reminiscence of the |A SC | 2 /T at zero temperature. 
However, a direct evaluation of the integrals yields, 


d£\d£ 2 


1 


7 T 2 y /1 + 441 + 4(41 + e 2 + 4l + £ 2 ) 


= 1 , 


d£\d£ 2 


41+4 + \/l + e 2 


l £ l| 


= - 1 . 


£ 2 


(37) 


Therefore, at leading order, Jn^ o (Q y ) = 0. This reproduces our previous result on the cancellation of leading order 
couplings between A^ and A sc . 

Now we repeat the above procedure for A c , the odd component of A]?- The bubble is given by 

d 2 k 

n fc 0 (Qy) = “ T > ( / — Vi Tr [^i (u m , k) ■ t 0 ■ g 2 (ujm, k ) • r 0 ] 


= -r£ 

m * 

= 2T E / 


1 ( 2 *)= 

d 2 k u 2 m - ei(fc)e 2 (fc) - |A SC | 2 
(27 rY^l + Elm^l + EKk)}- 


(38) 
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The only difference between Eq. (1551) and (1331) is the sign in front of |A SC | 2 . Subtracting the A sc -independent part 
and rewriting, we obtain 


'SnL(Q) = 


de\de 2 


|A a 


|A S 


8n 2 v x Vy [ EiE 2 {Ei + E 2 ) \Ei + E 2 |ei| + |e 2 

f°° deid£ 2 1 


8 V X Vy 


7T ~ \/l + £ i\/l + £ | (\/l + e i + \/l + e 2 ) 


+ 


deide 2 


7T 2 \ y /1 + e\ + y/1 + e\ ki| + ka| 


4' v x v v 


(39) 


where in the 


last line we have inserted the value of the two integrals in Eqs. (1571) in the bracket. 


B. The residual coupling between density component of charge order and superconductivity 

Since the leading order coupling between Ad and A sc vanishes, we need to go beyond leading order to compute their 
“residual” coupling, and we need to verify that this coupling is indeed repulsive. This can be done if one explicitly 
keeps the upper momentum cutoff A in evaluating /3 h and /3i 2 . However the analytical computation for the /3’s with 
a circular hard or soft momentum cutoff is more difficult than the computation of the bubbles in the normal state. 
Here we consider a case when a cutoff Ai, A 2 is imposed on e\ and e 2 respectively. We have evaluated the /3’s with 
more realistic cutoff’s numerically and have found that the result is qualitatively the same. 

We start from Eq. m- After converting the integration variable from (k X7 k y ) to (ei, e 2 ), and imposing a integration 
cutoff as described above, we have, 


Pn =2 TJ2 

m 

ih2 = -Tj2 


r A 


1 f A2 de\de 2 oj. 


m 


l 


-Ai J- a 2 8n 2 v x v y (w 2 * + el) 2 u4 + e\ ' 

Ai 


de\de 2 


1 


1 


I -Aj J- a 2 8tt 2 v x v v + ef + e\ 


2 ' 


We integrate over ei i2 first, 


(40) 


Pn = 


T 


27 J 2 V X Vy 


E 


■ tan 


^ 12 27 T 2 v x v y E 


■ tan 


-l 




_Ai 

| Urr, 

Ei 

I OJrr, 


Ai 


A? + W 2 , 
1 


■ tan 


■ tan 


-l 


UJr, 


a 2 


Wn 


Ao 


(41) 


The major contributions to both Matsubara summations come from the infrared, therefore we can expand the sum¬ 
mand as Laurent series of |w m |/Ai j2 . Up to next-leading order, we have 


Pn 


Pl 2 


T 


27 TV X Vy 

-T 

2-KVxVy 


y (-) 2 — 

m l Wrrl 

Efn 2 — 

^ V2/ \u m 


7T 


'm| 2A 2 |07 m | 

7T /I 


1 


|^m | 2 2|u7 m | \ Ai A; 


(42) 


where • • • stands for terms that are not divergent in l/|w m |. 

The first terms in both summands give the same results as shown in Eq. (1551) . which cancel out in (/3 n + Pi 2 ). 
Therefore, the bi-quadratic coupling between A^ and A sc is 


Pu + P12 = 


1 


log 


A 


47 tAi v x v y \T 


(43) 
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while the bi-quadratic coupling between A c and A sc , at leading order, is 


/5n — P 12 


1 

16 v x v y T' 


(44) 


Eq. (1431) is consistent with the experiment^, since it gives rise to a weak but finite suppression of the CDW correlation 
length upon entering the superconducting phase. 
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